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Abstract 

In the newly discovered iron-based superconductors, many experiments have demonstrated the existence of the ro- 
tational symmetry breaking nematic order, which has been a prevailing phenomenon in many correlated electronic 
systems. In this paper, we review nematic behaviors in iron-pnictides and the mechanism behind the development of 
the nematic order. We discuss evidences that support the spin-driven nematicity in iron-pnictides. Theories, results 
and predictions will be discussed based on this picture. We also briefly discuss the generalization of this theory to the 
nematicity in iron-chalcogenides. 
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1. Introduction 

Microscopically, strongly correlated systems are usu- 
ally described by extended Hubbard models with pa- 
rameters such as Hubbard interaction, Hunds coupling, 
doping, etc. Despite the relatively simple microscopic 
models, the infrared physics of strongly correlated sys- 
tems can be incredibly rich, various phases with com- 
pletely difl'erent low energy descriptions can appear in 
the same phase diagram tuned by a few microscopic tun- 
ing parameters. Directly deriving the precise low energy 
phase diagram from the microscopic extended Hubbard 
models is usually tremendously difficult except for some 
special limits of the model. However, the competition 
and interplay between low energy phases can still be 
reasonably well described using eff'ective theories that 
are based on some fundamental rules of physics. Many 
qualitative and semi-quantitative conclusions can be al- 
ready drawn based on these low energy eff'ective theo- 
ries without the detailed information of the microscopic 
physics of the system. 

Two decades after the discovery of the Cu-based su- 
perconductors (curpates), a completely new class of 
high temperature superconductors was found, where 
the Fe^^ ions play the most important rolef ll 0, [l, 
B B H [tI]- Iron-based superconductors can be di- 
vided into two basic categories: iron-pnictides and iron- 
chalcogenides. The category of iron-pnictides includes 
many diff'erent families of materials which are named as 



the 1 1 1 1-family (example: LaOFeAs) (EBItI], the 1 11- 
family (example: NaFeAs), the 122 family (example: 
BaFe2As2) Isl], and more complicated structures such as 
Sr4V206Fe2As2ll2]- Iron-chalcogenides include the 11 
family Fe/Se(Te)|[i3l and the 122 family-alkaH doped 
Ai_xFe2-ySe2 (A - K, Cs, Rb) llll IjJ, UJ, IjJ. The 
transition temperatures of these new superconductors 
reach as high as 56K||3) for iron-pnictides and 48K L15J 
for iron-chalcogenides. 

The discovery of Fe-based superconductors provides 
a precious opportunity to investigate possible universal- 
ities shared in all of high temperature superconductors. 
Fe-based superconductors share many common features 
with cuprates. Both of them are quasi-two dimensional 
systems with a similar doping-dependent phase diagram 
involving strong antiferromagnetism (AF)|16]. Super- 
conductivity is developed by introducing doping or ap- 
plying pressure while long-range magnetic orders are 
suppressed (for a review, see Ref. fPl]). Understanding 
the magnetism is thus pivotal to illustrate superconduct- 
ing mechanism in these materials. While the relation 
between magnetism and superconductivity has been a 
central research focus in the past several years, here we 
review another important aspect of magnetism: electron 
nematicity, in undoped or under-doped iron based su- 
perconductors. 

An electronic nematic liquid phase is characterized 
by rotational symmetry breaking. In the past decade. 
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this phase has been shown to appear almost universally 
in quasi-two dimensional electronic systems. The ne- 
maticity can arise from weakly correlated Fermi liquid 
systems dominated by kinetic energy to strongly corre- 
lated electron systems where electron-electron interac- 
tions are important. It has been observed throughout 
a variety of different electron systems including two- 
dimensional electron systems with high magnetic fields, 
strontium ruthenate materials and cuprates (for a review, 
see Ref. iH). 

In undoped or under-doped iron-pnictides, trans- 
port measurements reveal the existence of an intrinsic 
anisotropy of in-plane resistivity above magnetic Neel 
transition temperature 12111 . a strong evidence for 

the presence of an electronic nematic state. Similar 
anisotropy appears in a variety of experimental mea- 
surements, including local electronic structures mea- 
sured by scanning tunneling microscop y (STM) [221 , 



magnetic fluctuations by neutron scattering ||23l|24|], dy- 
namic conductivity by optical reflection measurements 
1I25I and electronic structures measured by angle re- 
solved photoemission spectroscopy (ARPES) 123] 
in de-twinned samples. 

However, it is still controversial whether the observed 
anisotropics indicate an indisputable electronic nematic 
state. The debate arises because of the interplay be- 
tween multi-degrees of freedom and the complexity of 
electronic structures in iron-pnictides. The lattice, spin, 
and orbital degrees of freedom all manifest themselves 
in the nematic state, which causes a debate on the origin 
of the rotational symmetry breaking of the tetragonal 
lattice. The magnetically ordered state in iron-pnictides 
is a collinear-AF (C-AF) state whose ordered wavevec- 
tors are (0, tt) or (tt, 0) with respect to the tetragonal iron 
lattice. The C-AF state not only breaks the S0(3) sym- 
metry of the spin space, but also breaks the rota- 
tional symmetry of the tetragonal lattice. Therefore, it 
is possible that an electronic nematic order which only 
breaks the latter symmetry can be induced by mag- 



netic fluctuations 1281 l22l]- Neutron scattering experi- 
ments have provided strong evidence showing that the 
magnetic fluctuations are indeed nematic above rA^| 24]. 
However, the rotational symmetry breaking can also 
be caused by other degrees of freedom. First, in iron- 
pnictides, the development of the C-AF order is always 
accompanied with a lattice structure distortion lUil Isill . 
A structural transition temperature Tg from the tetrago- 
nal lattice structure to an orthorhombic lattice structure 
is either at or above T^, namely, Tg > TnM]^^- Sec- 
ond, an onsite ferro-orbital order between dxz and dy^ 
orbitals also breaks the same rotational symmetry. 
Recently, such ferro-orbital order has been observed by 



ARPES E 12711 in pressure de-twinned samples. From 
standard symmetry argument, all the orders that break 
the same symmetry are allowed to be coupled directly 
with each other. Thus the existence of any one of them 
could lead to the presence of the others. Therefore, it is 
difficult to disentangle the role of the diff'erent degrees 
of freedom behind the rotational symmetry breaking. 

This paper is to review the progress based on the spin- 
driven nematic order in Fe-pnictides materials. Before 
we start our review, we would like to discuss accumu- 
lative evidence supporting that the spin degree of free- 
dom plays the dominant role in driving the rotational 
symmetry breaking. First, the lattice distortion observed 
throughout diff'erent families of iron-pncitides is rather 
small, much less than 1% of original lattice|[l3, Isill . 
Such a small lattice distortion can not cause a large 
change in electronic structures observed in the ne- 
matic state. Second, ferro-orbital ordering observed in 
ARPES ||26[ |27[ [32I1 is also small and a separated ferro- 
orbital ordering transition is never observed. More im- 
portantly, the change of electronic structure in nematic 
state is not dominated by the dxz and dy^ orbitals which 
can cause a rotational symmetry breaking. Instead, 
all the orbitals act together and are equally important in 
the process of band reconstruction observed in the ne- 
matic state| 331. Third, even if the ordered magnetic mo- 
ment in the C-AF state varies significantly in diff'erent 
families of iron-pnictides, the total magnetic moment 
obtained from inelastic neutron scattering is larger than 
over 2.0yUfi in underdoped samples[34, 24]. There is no 
question that the fluctuation of such a large moment can 
cause an electronic nematicity. Finally, as we will dis- 
cuss in the paper, predictions derived from spin-driven 
mechanism are in a good agreement with experimental 
resultsf29",'28]. 

Iron-pnictides are metallic instead of insulating. 
However, electron-electron correlation has been proved 
to be very important. Experimentally, it has clearly 
shown that the nematicity in iron-pnictides is not a 
Pomeranchuk type instability characterized by Fermi 
surface deformation. Instead, in iron-pnictides, the 
entire band structure is reconstructed in the nematic 
state I35I [33I1, which suggests that it is not driven by 
the states near the Fermi surface, rather by all electrons 
due to correlation efl'ects. These two facts lead to dif- 
ferent approaches to study spin-driven nematicity. On 
one hand, being a metallic state, it is not surprising to 
treat the C-AF state as a spin density wave(SDW) state 
and formulate the nematicity based on the two SDW 
order parameters specified by order wavevectors (0, tt) 
and (tt, 0) respectively! 37, 38]. On the other hand, due 
to the importance of electron-electron correlations as 
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Figure 1: (a), the lattice structure of Fe-pnictides materials at room 
temperature. The gray circles are Fe ions, the green and red circles 
are As ions above and below the Fe plane, respectively. The dashed 
square is the two-Fe unit cell, (b), the lattice structure after the Ising 
nematic order is developed: the thick lines represent the bonds be- 
tween parallelly aligned spins, but an Ising nematic order does not 
necessarily imply a nonzero long range spin order 1 16]. (c), the (0,;:) 
spin order at low temperature, (d), the sketchy phase diagram for lat- 
tice distortion and C-AF (or SDW) as a function of temperature and 
doping X, for some of the Fe-pnictides materials. The blue curve rep- 
resents the Ising transition, the red curve represents the C-AF transi- 
tion. 



well as the large magnetic moments observed in neu- 
tron scattering, we can take a local moment to for- 
mulate microscopic models based on short range ex- 
change interactions ||29[ |28|] . The nematic transition can 
be explicitly obtained and the nematicity can also be 
reflected in the magnetic exchange couplings. Fortu- 
nately, these two approaches result in an identical low 
energy Ginzburg-Laudau efl'ective action. 

In this review paper, we will first take a phenomeno- 
logical approach. We will discuss the field theory of 
the coupled nematic order and various aspects of the 
system at finite temperature, in particular, the nature 
of the phase transitions. Second, we will analyze mi- 
croscopic models based on eff'ective magnetic exchange 
couplings to obtain quantitative results and discuss the 
microscopic parameters that control the nematic transi- 
tions. Finally, we will discuss briefly nematic aspects of 
iron-chalcogenides, although at this stage it is still too 
early to write a complete review on these materials. 



2. field theory of nematicity 

2.7. Ginzburg -Landau field theory 

Based on the observation discussed in the previous 
section, we first seek a phenomenological field theory 
description of the (tt, 0), (0, n) C-AF (or SDW) and the 
nematic order. It was first proposed in Ref . list IS US] 
that the minimal spin model for the undoped parent Fe- 
pnictides material is either a5' = lor5'=2 spin model 
with nearest and next-nearest neighbor couplings Ji , J2 
that depend on the competition between the onsite Hub- 
bard interaction, crystal field splitting, and the Hunds 
rule: 



(1) 



<l,J> 



Unlike cuprates, in this system J2 is comparable with 
Ji because the As (or P) atoms are located at the pla- 
quette center of the square lattice of magnetic Fe ions, 
rather than on the links. There is also a much weaker in- 
terlay er coupling Jz, which is necessary to stabilize the 
spin order. It is well-known that when Ji < 2 J 2, the 
classical ground state manifold of the model in Eq.[T]is 

® S^, because the two sublattices of the square lat- 
tice will each form a Neel order (0i and 02 ), and the 
ground state energy is independent of the relative an- 
gle between these two Neel vectors. However, there 
is no generic symmetry of the system that protects the 
degeneracy within the classical ® ground state 
manifold (GSM), thus quantum and thermal fluctuations 
will both lift the degeneracy, leading to parallel or an- 
tiparallel alignment of the two sublattice Neel vectors 
|48, 49, 50], thus the generic GSM of the system is 

Z2, where the Z2 can be described by the Ising 
order parameter cr defined as follows: 



cr ~ 01 • 02. 



(2) 



We will defer the quantitative discussion of Eq. [T] 
to the next section, right now we will just discuss a 
Ginzburg-Landau theory for this lattice model. The lat- 
tice structure at high temperature has a perfect reflection 



symmetry 



xy 



y 



X, and under this trans- 



formation, cr always changes sign: cr -cr. Thus if 
(cr) 9^ 0, the reflection symmetry F^y has to break, i.e. 
the system symmetry must reduce from tetragonal to or- 
thorhombic. cr is precisely our nematic order parameter. 

Based on the symmetry of the system, the minimal 
GL theory describing the C-AF and nematic order is 



Fgl = (V;,cr)2 + r^cr^ + ^(V^0j2 
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Figure 2: The schematic phase diagram of Ginzburg-Landau mean 
field theory in Eq. [3] plotted against r = r^^ + r^, and Ar = - r^^. 
r is linear with temperature T, while Ar is tuned by anisotropy ratio 
JzlJin- When Ar is small, the interaction between 0i and $2 induces 
a strong first order transition, which corresponds to the undoped 122 
materials with more isotropic electron kinetics; when Ar is large, the 
transition is split into two transitions, with an Ising transition followed 
by an magnetic transition at lower temperature, and this is the case in 
the 1111 materials with quasi two dimensional physics. 



(3) 



Roughly speaking, r^- and are tuned by the transition 
temperature for nematic and C-AF order respectively, 
large r corresponds to the high temperature disordered 
phase. The cubic coupling u in Eq.[3]implies that when 
01 and 02 are ordered, the nematic order cr has to order; 
while the reverse is not necessarily true: there can be a 
phase with cr ordered while 0^ disordered. 

By simply minimizing the GL theory, we obtain the 
phase diagram in Fig. O When Ar = - r^- is large, 
meaning a has a much stronger tendency to order com- 
pared with 0^, the system goes through two second or- 
der transitions associated with ordering of cr and 0« re- 
spectively; while when Ar is small, the cubic coupling 
u in Eq.[3]merge the two transitions into one strong first 
order transition. 

What does Ar correspond to in the real system? It is 
reasonable to start with a quasi two dimensional struc- 
ture with interlay er spin coupling much weaker than 
inplane couplings /i ~ /2 ~ J- In Fe-pnicdides materi- 
als, the inplane coupling / is the energy scale that con- 
trols the ordering temperature of the Ising order param- 
eter cr. However, for purely two dimensional systems, 
the 0(3) C-AF is destroyed by thermal fluctuation at in- 
finitesimal nonzero temperature, thus there is only one 
Ising transition at finite temperature for pure two dimen- 
sional system. With small the transition temperature 
of magnetic order scales as Tj^ ~ J/(logJ/Jz) <^ T ising- 
Thus Ar is actually tuned by the eff'ective dimensional- 



ity of the system: When the system is more isotropic, 
Ar is smaller. The phase diagram of GL theory Eq.[3]is 
shown in Fig.|2l 

Experiments have shown that in the 1111 materials, at 
finite temperature there are two second order transitions 
associated with lattice distortion and magnetic transi- 
tion respectively ll39l] : while in undoped 122 materials 
there is only one strong first order transition where both 
phenomena occur simultaneously [4^ EH 112, H Hi . 
A simple Comparison between the experimental facts 
and our GL theory would imply that the 122 materials 
are more isotropic than 1111 materials. This is indeed 
true in real systems: First of all, the electron band struc- 
ture calculated from LDA shows a much weaker z di- 
rection dispersion compared with the 122 samples |51]; 
also the upper critical field Hd of 122 samples is much 
more isotropic [52]. This justifies treating the 1 1 1 1 ma- 
terials as a quasi two dimensional system, while treat- 
ing the 122 materials as a three dimensional one. When 
Jz and / are close enough, Ar is small, and the cou- 
pling between the nematic order and the C-AF will drive 
the transition first order by minimizing the free energy 
Eq.[3l 



2.2. finite temperature phase transitions 

Now let us focus on the regime with large Ar, i.e. 
when there are two separate second order transitions 
at finite temperature. In principle, the transition of the 
Ising order parameter cr should belong to the 3D Ising 
Wilson-Fisher (WF) universality class, while the transi- 
tion of the C-AF order parameter should belong to the 
3d 0(3) WF universality class. Thus at both transitions 
the specific heat should have a singularity peak. How- 
ever, the specific heat measurement on BaFe2-xCoxAs2 
reveals two close but separate transitions at finite tem- 
perature, with a sharp peak at the C-AF transition, and 
a discontinuity at the lattice distortion transition ll53ll . A 
discontinuity of specific heat is a signature of mean field 
transition, in contrast to the sharp peak of Wilson-Fisher 
fixed point in 3 dimensional space. 

Mean field transition can only be accurate for sys- 
tems with spatial dimension equal to or higher than 4. 
Thus the specific heat measurement in Ref . ll53n implies 
that eff'ectively the dimension of the nematic transition 
is enhanced. In this section we will demonstrate that by 
coupling to the lattice strain field fluctuation, the efl'ec- 
tive spatial dimension Deff for the nematic transition is 
enhanced to Deff = 5, while the dimension for C-AF 
transition remains Deff = 3. 

Let us first investigate the C-AF transition. The C-AF 
transition at finite temperature should belong to the 3D 
0(3) transition if the lattice elasticity is ignored. The 
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0(3) order parameter couples to the lattice strain field 
with a quadratic term 15411 : 



^(dxUx + dyUy + A^d^Uz). 



(4) 



After integrating out the displacement vector w, this cou- 
pling generates a singular long range interaction be- 
tween |0p in the real space: 



- r'|3 



(5) 



/ is a dimensionless function which depends on the di- 
rection of r - r^. The scaling dimension of g is A[g] = 
a = 2/y - 3, and v is the standard exponent at the 3D 
0(3) transition, which is greater than 2/3 according to 
various types of numerical computations |55]. There- 
fore this long range interaction is irrelevant at the 3D 
0(3) transition, and by coupling to the strain field of 
the lattice, the C-AF transition is unaff'ected. However, 
if the C-AF has an Ising uniaxial anisotropy, the C-AF 
transition becomes a 3D Ising transition with y < 2/3, 
and the strain field would induce a relevant long range 
interaction [ 54], which leads to a run-away flow from 
the 3D Ising fixed point. 

However, since the symmetry of the nematic order 
parameter cr is the same as the shear strain of the lattice, 
the strain tensor will couple to the Ising field cr as |[3Q|] 



= Acr(dxUy -\- dyUx) -\- • • 



(6) 



it is the displacement vector, the ellipses are all the elas- 
tic modulus terms. Notice that we have rotated the co- 
ordinates by 45 degree, since the true unit cell of the 
system is a two Iron unit cell (Fig.O. After integrating 
out the displacement vector it, the eff'ective free energy 
of cr gains a new singular term at small momentum: 



(7) 



/ is a function of spherical coordinates 6 and defined 
as (kx,ky,kz) = ^(cos(^) cos(0), cos(^) sin(0), sin(^)), 
but / is independent of the magnitude of momentum k. 
By tuning the uniform susceptibility r, at some spherical 
angle the minima of / reduces to zero, we will call these 
minima as nodal points. These nodal points are isolated 
from each other on the two dimension unit sphere la- 
belled by 6, 0, and are distributed on the unit sphere 
(6, 0) according to the lattice symmetry. 

The symmetry allows the minima to locate at the 
north and south pole, i.e. 6 = and tt. Expanded at 
these two minima, the free energy reads 



The naive power counting shows that eff'ectively the 
spatial dimension of this field theory Eq. [8] is De/f = 
5. The quartic coupling cr^ is irrelevant based on the 
straightforward power counting. Therefore the strain 
tensor fluctuation efl'ectively increases the dimension by 
two, which drives the transition a mean field transition. 
This result is completely consistent with the specific 
heat measurement on BaFe2-xCoxAs2. 

Using the Cartesian coordinate, the free energy Eq. [8] 
becomes: 



F= \ dkxdkydU^^^+kl + r)\(Tk\' + (9) 



Based on this equation, the eff'ective dimensions for mo- 
menta kx, ky and kz are 



Afc] = A[/:J 



2A[/:J 



(10) 



(8) 



Thus the total dimension is still 5, considering /S\kx\ = 
A[ky] = 2A[^^] = 2. All the other momentum dependent 
terms in this free energy are irrelevant. 

Based on the discussions above, we conclude that the 
lattice elasticity fluctuation will strongly modify the na- 
ture of the nematic transition, but will leave the C-AF 
transition unafl'ected. In Ref. (s^, the authors proposed 
that when the two transitions are very close, the lattice 
fluctuation has the potential of driving the C-AF transi- 
tion first order, because when these two transitions are 
very close the eff'ective quartic term of the magnetic or- 
der parameter can become negative. This eff'ect was also 
observed in experiments on lightly electron doped 122 
materials in Ref. |20]. 

2.3. Hertz-Millis theory 

In this subsection we will discuss the quantum critical 
points associated with nematic and C-AF orders i3Qll. 
A Hertz-Millis theory for quantum nematic order was 
studied in Ref. |56]. However, in our system the ne- 
matic order is strongly correlated with the spin interac- 
tion, thus the Hertz-Millis theory will be very difl'erent. 
The LDA calculation and photoemission both conclude 
that in the Brillouin zone of most Fe-pnictides materi- 
als, there are two hole pockets at the F point (0, 0), and 
two electron pockets located at the M points ((0, n) and 
(tt, 0)). Since the hole and electron pockets have difl'er- 
ent shapes, under translation of 2 = (tt, 0) in the mo- 
mentum space, the hole pocket will intersect with the 
electron pocket. This intersection leads to overdamping 
of the SDW order parameters. The decay rate can be 
calculated using Fermi's Golden rule: 



OJ 



~ Co— — • 

\Vh X Ve\ 

Vh and Ve are the fermi velocity at the points on the 
hole and electron pockets which are connected by wave 
vector (tt, 0). The standard Hertz-Millis |57] formalism 
leads to a coupled z = 2 theory in the Euclidean mo- 
mentum space with Lagrangian 
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Figure 3: The solution of the RG equation Eq. [161 AH three quar- 
tic perturbations decrease first, then increase and finally become non- 
^/2_ 2\ ^ perturbative. The run-away flow is weaker than marginally relevant 

7yI\Qx ^y) r2 'perturbations. 



L = A(|0i|4 + |02l')-a(0i.02)' + C|0ip|02p. (12) 

The parameter r can be tuned by pressure and doping. 
The Ising symmetry of cr = 0i • 02 on this system corre- 
sponds to transformation 



perturbati( 

dB 



01, 02 ^ -02, cr -cr. 



(13) 



This Ising symmetry excludes the term 0i -02 in the 
Lagrangian, while the mixing term y(pi(ql - q^) • 02 is 
allowed. 

We can diagonalize the quadratic part of this La- 
grangian by defining 0a = (0i + 02)/V2 and 0^ = 
(0i-02)/V2: 



H ' (\0J\ + (1 - ^)qi + (1 + ^)q'y + r)0A 



+ $B • (1^1 + (1 + \)ql + (1 - \)ql + r)$B + L\ 

L = A(|0A|4 + |05|4) + 5(0A-05)' + C|0Ah05a4) 
After the redefinition, the Ising transformation becomes 



0fi, 0fi -0A 



cr -cr. 



(15) 



Naively all three quartic terms A, B and C are marginal 
perturbations on the z = 2 mean field theory, a cou- 
pled renormalization group (RG) equation is required to 
determine the ultimate fate of these terms. Notice that 
the anisotropy of the dispersion of 0a and 0^ cannot be 
eliminated by redefining space and time, therefore the 
number y will enter the RG equation as a coefficient. 
The final coupled RG equation at the quadratic order 
for A, 5 and C reads: 



dA 
d]ni 



-22A^ 



1 



-B' 



-C^-BQ 



d\nl 
dC 
dhd 



-5uB^ - MB - SuBC, 



-uB^ - 4AB - 20AC - 4uC\ 



(16) 



w is a smooth function of y, which decreases smoothly 
from w = 1 in the isotropic limit with y = 0, to u = 
in the anisotropic limit with y = 2. The self-energy 
correction of 0^ from the quartic terms will lead to the 
flow of the anisotropy ratio y under RG, but this flow is 
at even higher order of the coupling constants, thus in 
our current discussion u is just taken a constant. 

The typical solution of the RG equation Eq. [16] is 
plotted in Fig.[3l for the most natural choice of the initial 
values of parameters with A » a, A » C in Eq. \T2\i.e. 
the coupling between 0i and 02 is weak. One can see 
that the three parameters A, B and C all have run-away 
flows and eventually become nonperturbative. This run- 
away flow likely drives the transition weakly first order. 
However, the three coefficients will first decrease and 
then increase under RG flow. This behavior implies that 
this run-away flow is extremely weak, or more precisely 
even weaker than marginally relevant perturbations, be- 
cause marginally relevant operators will still monotoni- 
cally increase under RG flow, although increases slowly. 
Therefore in order to see this run-away flow, the correla- 
tion length has to be extremely long i.e. the system has 
to be very close to the transition, so the transition re- 
mains one single second order mean field transition for 
a very large length and energy range. 

The interlay er coupling of 0a and 0^ has so far been 
ignored, this is also a relevant perturbation at the z = 
2 mean field fixed point. The z direction tunnelling is 
written as Jz^a,z ' 0«,z+i, which has scaling dimension 2 
at the z = 2, J = 2 mean field fixed point, and it becomes 
nonperturbative when 



(17) 



This equation implies that if the tuning parameter r is 
in the small window r < Jz/Jin, the transition crossover 
back to a z = 2, J = 3 transition, where all the quartic 
couplings are irrelevant. Since in the two dimensional 
theory these quartic terms are only very weakly relevant, 
in the end the interlayer coupling may win the race 
of the RG flow, and this transition becomes one stable 
mean field second order transition. 



3. Microscopic theory 

3.1. Model 

As we have mentioned earlier, a simple microscopic 
model to capture the C-AF state is the Ji -J2-Jz Heisen- 
berg model. The interlayer coupling must be included 
in order to address finite temperature magnetic transi- 
tion. Moreover, as we will show later, the layer coupling 
Jz provides a control of the diff'erence between nematic 
and magnetic transition temperatures. The Hamiltonian 
for such a model can be written as, 

n {ij) m) 

-\-Jz ^ S/,„ • S/,„+i], (18) 

n,i 

where n is layer index. The spin nematic order which 
breaks the rotational symmetry can be defined as 



(19) 



where x and y are unit vectors in a two dimensional 
tetragonal plane. 

As we have mentioned in the introduction that both 
lattice distortion and ferro-orbital ordering which break 
the same symmetry can be coupled to the spin ne- 
matic order. To be more specific, we consider the NN 
coupling /i(|r/ - r^l) as a function of the lattice distance 
between two sites. If there is a small orthorhombic lat- 
tice distortion Sr, the NN coupling becomes 



/i(ro + Sr/2)SiSi^, + /i(ro - ^r/2)S,-S, 



i+y 



(20) 



where ro is the lattice constant of the original tetrago- 
nal lattice. We can expand the above term and obtain a 
coupling between the lattice distortion and the nematic 
order, Asi(ro)Srcri, where Asi(r) = There is also 

coupling between the spin nematic order and an onsite 
orbital order between dxz and dyz orbitals. The coupling 
can be written as 



(n^zii) - nyz{i))cr, 



(21) 



where nxz(yz)(i) is electron density at the corresponding 
orbitals. If we include all these couplings and consider 
that the magnetism is the dominant player in the sys- 
tem, we can integrate out the above lattice and orbital 
fluctuations to obtain an efl'ective pure magnetic term to 
account for the couplings. It is easy to show that such an 
integration leads to a biqradratic spin coupling limited 
to the NN bonds, namely. 



(22) 



<iJ>NN,n 



By adding such a term into the Hamiltonian in EqfTSl 
which will be called as the Ji - J2 - Jz - K model. 
The C-AF state is the ground state of the model when 
Ji < 2/2 and ^ > 0. This extension has been recently 
considered in 15^ ^Q] . 



3.2. Large-S limit 

We can exactly solve above model in the large S- 
limit. Without the K term, in the case of Ji < 2 J 2, 
the ground state is infinitely degenerate, as the relative 
angle between the spin directions on two sublattices is 
entirely arbitrary. If we consider quantum fluctuation, 
such a degeneracy will be broken. To show this, we can 
take Holstein-Primakofl' transform for spin operators: 



^r = V2^«,Vi-^,5-=V2^«j^i-^ 



S] = S - a] at, 

where a/'s are annihilation operators of spin wave 
bosons satisfying [auaj] = [aj,ap = and [a/,ap = 
Sij. We consider the following spin configuration in 
two sublattices: if a site (ix, iy) belongs to the first sub- 
lattice, its spin moment is Sij = 5'(0, 0, (-1)^0; if it 
belongs to the second sublattice, its spin moment is 
= 5'(-l)^'^(sin^, 0, cos^). Taking linear spin wave 
approximation in the Holstein-Primakoff' transform, we 
have 

S u = V25 ( ai )«; , 



(23) 



(24) 



V25"( 



l+(-l)'- 



1 -(-!)'« 



)ai. 



= (- ly '( -yj - (ai +01) cos + {S -ajai) sin 6) 



(25) 



2,1 



(-!)''( -J— (a; + al) cos 6 + {S -aloi) sin 6) 
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-(ai-aj). 



(26) 



'2,i 



(-ly^C- + (^J) sin 6 + (S - a] ad cos 6). 



Substituting these expressions into Eq.([T8l), and taking 
Fourier transformation, we obtain 

^Lsw Y^YAkaluk + ^{aua-k + h.c.)\, (27) 



where 



200 



> 



B'lOO 




1.0 1.5 
(1,K,1)(r.l.u.) 



Ak = 4JiS + 4J2S + 2JiS [(1 + cos 0) cos(k^) 

+ (1 - cos 0) cos(ky)] + 4/„ (28) 

Bk = -{2/i^ [(1 - cos 6) cos(k^) 

+ (1 + cos ^) C0S(^3;)] + 2/2 5" (COS(^;c + ky) 

+ cos(^;, - ^3;)) + 4J,cos(k,)}. (29) 

Taking the BogHubov transform for bosons, we ob- 
tain 



H 



where 



LSW 



(30) 



(31) 



and jk = coshi/zkak + sin\n//ka[ 



cosh i//kCil 



sinh i//ka-k, and tanh(2j/ry^) = At zero temperature 
for small the quantum fluctuation contributes a zero- 
point energy for bosons 



7Q 



IE 



cok = -^^rr^ cos^ + J20{{^)\ (32) 



where jq = 0.13. The ground state, therefore, is either 
^ = or ^ = TT, which means that spins in two sublat- 
tice are collinearly aligned, namely, the C-AF state. The 
collinear configuration is an eff'ect of spin fluctuation. 
This eff'ect is called "order by disorder" 14811491 |50ll. 

The energy in Eq(32lcan be taken into account by the 
bi-quadratic K term as well. The quantum fluctuation 

contribution to the value of K is given by We 
can thus view the K term in the Ji - J2 - Jz - K model 
as an efl'ective term that includes contributions from all 
difl'erent couplings. In the C-AF state, the K term can 
be decoupled as 



Hq = -2 ^ Kcrij^n(Si,n ' ^ j,n) + COHS. 



(33) 



<iJ>NN,n 



where (T^„ =< S/,„ -S^^ >= ±cr, is the expectation value 
of the spin- spin interaction in the C-AF state. The sign 



Figure 4: The spin wave dispersion along (IT, k, IT) direction reported 
in|23]: the soHd line is fitted when Jia » J\b and the dashed line is 
fitted if J\a ~ J\b- 



is chosen to be positive if the spins are FM along one 
axis (b-axis) and negative if they are AF along the other 
axis (a-axis) in the C-AF state. Combining this term 
with J I term in the original Hamiltonain, the efl'ective 
nearest neighbour exchange couplings in the C-AF State 
become 



Jia = Ji-^ 2Kcr, Jiij = Ji - 2Kcr. 



(34) 



Thus, the spin wave excitations in the C-AF state are 
efl'ectively described by a Jia - Jib - J2 - Jz model. 
Experimentally, a large difl'erence between J] a and 



Jih is indeed observed in the C-AF state|6ll S S). 
The spin wave energy at wavevector (tt, tt, tt) is very sen- 
sitive to the difl'erence between Jia and Ji^. The en- 
ergy is proportional to ^JJla - Jib- Experimentally, The 
gap coin, 71, 7i) measured in the 122 family has been very 
large as shown in fig HI namely, the diff'erence of Jia 
and Jib is very large. Such a large diff'erence suggests a 
strong nematic order and a sizable K value. 

3.3. large -N limit 

The model can be analytically solved in the large-N 
limit to obtain phase diagram and transition tempera- 
ture. We take the continuum limit following the deriva- 
tion m The low energy eff'ective Hamiltonian 
reads 

Heff = I (fr^[J2\Vnc,nf - Jz^a,n ' n„+i,^] 

^ n,a 

-^^[ni,„ • n2,nf + y X ' (35) 

n n 

where // = JiS^, K = KS^/2 and ni(2),„ are AF orders 
in the two sublattices mentioned earlier respectively and 
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Figure 5: LeftiT^ and T^^ as the function of for J2 = 2Ji and N = 3. 
Right: ^'xn^ function of for Ji = O.5/2. 



the nematic order in the continuum Hmit becomes cr = 
• n2,«. We can solve the effective model in the large- 
N limit analytically. In general, there are two transition 
temperatures. The nematic transition temperature Ts is 
analytically given by 



47r/2 

at7 



In- 



JiKNTs) 



(36) 



^NTs'' ' 47TJ2 

and the C-AF transition temperatures, is determined 

by, 



— = — — In 
2K 87r/2 

CTc 1 

2I 



1 . /2 

— — In — . 

47r/2 /z 



(37) 



T]^ is always lower than Ts in the solutions of above 
equations. Moreover, = if = 0, and is finite 



Figure 6: Plot of structural phase transition temperature Ts, and mag- 
netic transition temperature Tj^ associated to the AF order of Fe ions, 
versus the c-axis for RFeAsO (R = La, Sm, Gd, Tb). Inset: plot of 
{Ts - Tn)/Tn versus Jz/Ji- The open circles denote the experimen- 
tal data of (Ts - T^)IT^, and the solid squares denote the theoretical 
values I63n . 



Otherwise. For small one can approximately have 
Tn ~ -^2/111(4). In fig 13 we plot the two transition tem- 
peratures and their diff'erence as the function of JzlJi- 
The diff'erence decreases as increases. This qualita- 
tive prediction is consistent with the GL field theory ap- 
proach discussed in the previous section, and it is also 
confirmed in the 1 1 1 1 family! 63] as shown in figlU 

This calculation shows that as soon as is much 
larger than 10"^/^, the two transitions are close to each 
other so that they become practically inseparable in ex- 
periments. This result is consistent with the fact there 
is only one phase transition in the parent compounds of 
the 122 families where determined by spin wave ex- 
citations is rather large: in CaFe2Se2, Jz is almost one 



third of the in-plane coupling /2|23]. In BaFe2As2, Jz is 
about 0.015/2124]. The prediction is further confirmed 
in under-doped 122 materials. Upon electron doping 
in BaFe2As2, for example, B aFe? _ v Niv As? |[6411 . the spin 
excitations become much less three dimensional. The 
eff'ective Jz is drastically reduced and two phase transi- 
tions were observed. 



4. Nematic order in Iron-Chalcogenides 

Iron-chalcogenides include two families, the 11 
family (FeTe/Se) and the 122 family (Ai_xFe2-ySe2 
(A = K, Rb, Cs). While the former shares similar band 
structures with iron-pnictides, the latter exhibits several 
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distinct characters that are noticeably absent in other 
iron-based superconductors, including the absence of 
hole pockets at T point of Brillouin zone in their su- 
perconducting (SC) phases |[65[ |66[ l67ll, AF ordered in- 
sulating phasesi68l[6l,[7^] with very high Neel transi- 
tion temperatures in their parental compoundsl 13] and 
intrinsic vacancy ordering. These intriguing properties 
have stimulated many exciting research activities. 

While it is too early to call a reasonable review 
for iron-chalcogenides, the magnetic orders in the par- 
ent compounds of these two families have been mea- 
sured. Unlike iron-pnictides whose different families 
share a common C-AF state, different families of iron- 
chalcogenides exhibit different magnetically ordered 
states. The FeTe compounds of the 1 1 family has a bi- 
collinear AF(B-CAF) ground state as shown in fig Ha). 
The B-CAF state is accompanied by a monoclinic lat- 
tice distortion. The distortion can produce 1.3% change 
of the tetragonal lattice constant, a value much larger 
than the orthorhombic distortion in iron-pnictides|71]. 
The large lattice distortion indicates a strong coupling 
between spin and lattice and is consistent with the fact 
there is only one strong first-order type transition in 
FeTe, namely, the B-CAF and monoclinic lattice distor- 
tion take place at the same transition temperature which 
is about The ordered magnetic moment in 

FeTe about 23jib is also significantly larger than the 
ones in iron-pnictides. 

For the 122 family, Ri_xFe2-3;Se2, one magnetically 
ordered state with iron vacancy order, called block- AFy 
has been identified when x = 0.8 and y = 0.4 as shown 
in fig|3b). Iron vacancies form sl ^^5 x ^^5 pattern and 
four iron spins group together to form a super- site with 
ferromagnetic(FM) alignment. In the lattice formed by 
the super-sites, it is AF ordered. The ordered moment 
per iron site in the B-AFy state reaches and the 

transition temperature is about 500k||68l[63]. The lattice 
distortion in the B-AFV state can reach 2% of the orig- 
inal lattice constant. More importantly, the B-AFy state 
is insulating with a gap larger than 0.3^v. It is the first 
insulating state observed in iron-based superconductors. 

Recently, it has been shown that the effective mag- 
netic exchange models that describe different ordered 
states can be unified into sl Ji - J2 - J3 - K - 
modelH H [60], where J3 is the next next nearest 
neighbour magnetic exchange coupling. Compared to 
the Ji - J2 - K - Jz model for iron-pnictides, there 
are two major differences: Ji is strong ferromagnetic 
(FM) in iron chalcogenides and J3 is significant and AF. 
Without vacancy, the model has two degenerate classi- 
cal ground states: the B-CAF or bloch-AF (B-AF) state 
which is shown in fig He). An effective theory based 



on this model has been discussed in [75]. The complete 
field description of the Ji - J2 - J3 - Jc - K model near 
the B-CAF or B-AF phase is given by Ht = Hb + Hk 
where 

Hb= f J^r^ [i/3|V4a(r)P - jJnA^) • 0«+i,a(r) 

-g' W • 43(r)]' + ' 44(r)]\ (38) 



and 



Hk = -K r J^r ^{|0,,i(r) . 0„,2(r)]' (39) 

-F[0„,i(r)-^„,4(r)f 
+ [0.,2(r) • 0„,3(r)]' + [0.,3(r) ' 0«,4(r)]'} (40) 



where $n,a with a = 1 , 2, 3, 4 are four Neel order param- 
eters defined in the four sublattices of the tretragonal 
lattice specified by the J3 coupling and g' ~ O.lSJ^/Js. 
In this model, there are four independent Ising or- 
ders CTi = g'$n,lOO • 0n,3(r)^ ^2 = ^'0«,2(r) ' 0«,4(r), 

cr+ = K{$n,i(x) + 0n,3(j)) • (0n,2(r) + 0«,4(r)) and cr_ = 
K{$n,i(x) - 0„,3(r)) • (0„,2(r) + $nA(x))' 0-+ and o-_ de- 
scribe the Ising orders associated with the B-AF and B- 
CAF phases respectively. The Ising order (Ti breaks ro- 
tational symmetry and represents a nematic phase. The 
analysis of this effective field theory has been carried 
out inlTSllfiS]- This effective theory also suggests that 
the lattice- spin coupling is extremely important in de- 
termining the true magnetic order in iron-chalcogenides 
due to intrinsic magnetic frustration. 

Since there are few experiments that have been im- 
plemented to probe the nematic nature of the B-CAF 
state explicitly, we will not extend our discussion further 
here. However, it is interesting to point out that since the 
122 family of iron-chalcogenides may approach a mag- 
netically ordered insulating phase, the materials may 
carry similar nematic phases as curpates upon doping. 
The charge degree of freedom hence may manifest it- 
self and play an important role in driving nematicity in 
these new materials. 



5. Summary and Extensions 

In this article we reviewed the study of the electronic 
nematic order and the C-AF order that have been univer- 
sally observed in iron-pnictides materials. Based on the 
experimental observation of the close relation between 
these two order parameters, we propose the theory of 
spin interaction driven nematic order in these materials. 
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of the lattice will significantly enrich the physics. 

The 122 family of iron-chalcogenides has opened a 
new chapter in the field. It provides a fresh ground to 
test the concept of nematism. Since electron-electron 
correlation is generally believed to be stronger in iron- 
chalcogenides than iron-pnictides, exciting results re- 
lated to nematism in the future can be expected. 

Acknowledge: We thank S. Kivelson, S. Sachdev, C. 
Fang, Pengcheng Dai, Donglei Feng and X.H. Chen for 
many valuable discussions. 
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Figure 7: Sketch of spin configurations of (a) bi-coUinear AF (B- 
CAF) state observed in FeTe, (b) block AF with vacancy order (B- 
AFv) observed in K0.8Fe1.6Se2, and (c) block AF order (B-AF) order. 



In this article we have thoroughly studied this problem 
with both low energy field theories such as Ginzburg- 
Landau theory and Hertz-Millis theory, and also micro- 
scopic theories in the large-^" and large-A/^ limit. We 
also briefly discussed the generalization of this picture 
to iron-chalcogenides materials. 

A proof of the nematic order induced by spin interac- 
tions can have strong impact on the field, it would auto- 
matically suggest that the superconductivity mechanism 
is most likely to be magnetically driven. While both 
experiments and theories reviewed here have provided 
convincing evidence, we agree that the debate over the 
origin of the rotational symmetry breaking has not been 
fully settled at this stage. 

Although we have a detailed theory of the interplay of 
electron nematicity and spin order, the relation between 
these order parameters and the high temperature super- 
conductor is still unclear. Strong competition between 
spin orders and superconductor phase has been observed 
experimentally in iron-pnictides, meaning that a com- 
plete description of low energy physics should also in- 
volve the superconducting phase. This is one of the di- 
rections that can be pursued in the future. 

In section II, we have seen that the lattice strain field 
fluctuation strongly afl'ects the finite temperature phase 
transition of the nematic order. However, in section II 
the discussion of the lattice strain field was purely clas- 
sical and phenomenological. At zero temperature, the 
nematic order will conceivably strongly couple to the 
phonon modes of the lattice, thus the phonon fluctuation 
will most likely also strongly afl'ect the quantum criti- 
cal behavior, as long as the spin wave excitation at zero 
temperature. A complete theory including the phonon 
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